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Abstract I study the sunspot area fluctuations over the epoch of 12 solar cycles 
(12-23). Lately, I found three significant quasi-periodicities at 10, 17 and 23 
solar rotations, but two longer periods could be treated as subharmonics of the 
10-rotation period. Thus, I search this period during the low- and the high- 
activity periods of each solar cycles. Because of the N-S asymmetry I consider 
^^' each solar hemisphere separately. The skewness of each fluctuation probability 

distribution suggests that the positive and the negative fluctuations could be are 
examined separately. To avoid the problem when a few strong fluctuations could 
create an auto-correlation or a wavelet peak, I also analyse the transformations 
^S| ' of fluctuations for which the amplitudes at the high- and the low-activity periods 

are almost the same. The auto-correlation and the wavelet analyses show that 
the 10-rotation period is mainly detected during the high-activity periods, but 
it also exists during a few low-activity periods. 
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1. Introduction 



^ ' In the last decades, the mid-term quasi-pcriodicitics of many solar activity trac- 

^D , crs have been discussed. IWolffI (|1983p reported the about 12 rotations periodicity 

in the monthly Wolf number variations from 1749 to 1979. This periodicity 
was also found in many solar activity indicators. |Delarche, Laclare and Sadsoud] 
(|1987p detected it in the solar diameter measurements during cycle 21. 
ILean and Bruecknerl (|1989p found it in the power spectrum of the sunspot block- 



^— N ' ing function, 10.7 cm radio flux, sunspot number and plage index daily data dur- 

ing the three solar cycles (Nos. 19-21). [Pap, Tobiska and Bouwer|(ll990p showed 
that an 8-11 month period existed in the total and UV irradiances (1980 and 
1982-1988 respectively), 10.7 radio flux (1947-1989), the Ca K plage index (1970- 

^^ ■ 1987), the sunspot blocking function (1874-1982) and the Mg ii core- wing ratio 

H I (1978-19861. lAkioka et al.\ (|1987p detected it in areas and numbers of sunspot 

groups from 1969 to 1986. lOliver et al.l (|1992p studied daily sunspot numbers 
from 6 to 21 solar cycles. They adopted two different power spectrum methods 
and found a periodicity of about 323 days in three solar cycles only. They also 



Astronomical Institute, University of Wroclaw, Poland, 
email: |getko@astro.uni.wroc.pl| 



SOLA: klu_02.tex; 24 May 2010; 0:34; p. 1 



R. Getko 

found a significant periodicity at 350 days (13 rotations) during cycles 12-21 to- 
gether and in some individual cycles. [Getkol (|2006p found statistically significant 
quasi-periodicity of about 9 months (10 rotations) in both high- and low-activity 
periods for the monthly Wolf number fluctuations from 1 to 22 solar cycles and 
for the group sunspot number fluctuations from 5 to 18 solar cycles. 

Two longer periods at 17 rotations and at 23 rotations were found in many 
solar activity parameters. More up to date review is by [Obridko and Shelting| 
(|2007p . Latelv. [Getkol (|2009p showed that they could be treated as subharmonics 
of the 10-rotation period. This facts could explain a wide range of periodicities 
in various solar indices at all levels from the tachocline to the Earth. 

In this paper, I search the 10-rotation period in the sunspot areas during 
the low- and the high-activity periods for cycles 12-23. The N-S asymmetry 
of solar activity (jVizoso and Ballesterl I1990p suggests that both hemispheres 
should be considered separately. Because each empirical probability distribution 
of fluctuations is asymmetrical, the positive and the negative fluctuations are 
considered separately. To avoid the problem when a few strong fluctuations could 
create an auto-correlation or a wavelet peak, I transform each fluctuation time 
series into a new time series which has a constant variance (jGetkoi r2006p . The 
amplitudes of each new time series at the high- and the low-activity periods are 
almost the same. The auto-correlation and the wavelet analyses of the original 
fluctuations and their transformations are used to find the mid-term periodicities 
from high- and low-activity periods in each solar cycle. 



2. The Stationary Version of Sunspot Area Fluctuations 

I consider the daily sunspot areas for the northern hemisphere (i*"), and the 
southern hemisphere (i^f) for solar cycles 12-23 available at the 
National Geophysical Data Center (NGDC) 
(http://solarscience.msfc.nasa.gov/greenwch/). For the i-th Carrington rotation 



I evaluate the mean area for the northern hemisphere (5f ): 5"," = -^ J2i=i ^l i 
where K is the number of days for the i-th rotation. I define the fluctuation (F") 
of the mean area (S^) from the smoothed mean area: F^ = 5" — S"" for i = 
1, . . . , TV, where ^=13 EjtLe ^l'- Each of {Pp} and {F/} contains N = 1706 
elements. The beginning and the end of each cycles is defined by the minima of 
the smoothed monthly Wolf numbers. Each time series has the zero mean value 
and the time-dependent variance. To stabilize the variance I apply the procedure 
which is shown by IGetkol (|2006p in Appendix. For the i-th rotation the value 
X" of the new time series {X"}^^ is a function of two parameters k and u 
(see equation (A2) in the Appendix). I estimate the values k and log A from 
equation (AID). The best solution is for m = 13 (rotations) and k = 0.68 ± 0.02. 
For {i^"} it is for u = 13 (rotations) and k = 0.83 ± 0.02. Both the coeflicients 
k are statistically significant (the F-statistic is used). 

To verify the stationarity conditions of {A"} and {A/} the means {X = 

^ ESll'tl ^^■)' ^t^^'l^^^ deviations (W= ^ iJ:^^l%iXj -Xf ), and 
auto-correlation functions (c^) are computed. For each time scries about 70% 
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Figure 1. (a) Auto-correlation functions (cr) of {X"}, {X"°} and {X^*} (for all the data 
combined). The dashed horizontal lines represent two standard errors for each Cr functions, 
(b) Same as for (a), but for {X|}, {Xf°} and {Xf*}. 



of the values of X for u = 13 and zi = 7, 20, . . . , 1698 belong to the interval 
(X - a,X + a), where X^ « -0.15 and (t" « 0.7, and A^ « -0.08, and 
a'^ Ri 0.3. All sample standard deviations for each time series belong to the 
interval (ct— 2ct, W+2a), where cr'^ w 4.4 and a" w 0.4 and cr® w 1.8 and 
P « 0.2. Figs. HJi andlB) show the c^ functions of {X,}, {X°} and {X*} for 
all 12 cycles for the northern and the southern hemispheres respectively. Here 
{X°} and {X*} are computed after dividing {Xi}^-^ into two parts: those from 
the periods of high and low sunspot activities, respectively. To evaluate these 
periods the smoothed monthly Wolf numbers {Rm}m=i ^'^s used. All m for 
which the values Rm are less (greater) than the mean value of {Rm} define 
low (high) activity periods. Two standard errors for each Cr are plotted as 
the dashed horizontal lines. Because each of three time series has a different 
length M, the bounds for the shortest time series are plotted (±0.068). The 
differences between the c^ functions of {X"°} and {X"*} are greater than 2(j 
for T = 11, 13, 16, 17, 20, 25, 33, 34 rotations. For the southern hemisphere both 
the Cr functions significantly differ at r = 7, 10, 12, 15, 19, 23, 30, 34 rotations. 
Although the means and the standard deviations of {ATi}, {X°} and {X*} are 
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Figure 2. (a) Top: Auto-corrclation function (cr) of {i^"} for cycle 18. Bottom Same as for 
tile upper curve, but for {X"}. Tlie dotted lines represent two standard errors of each Cr 
function, (b) The same as for (a), but for {i^" } and {X" }. (c) The same as for (a), but 
for {Ff-} and {Xf"}. 



aproximately constant the Cr functions are different. Thus, {Xf} and {Xf} are 
not stationary. 

Because each fluctuation probability distribution has positive skew (|Getkol 
I2009p . the positive and the negative fluctuations arc considered separately. For 
{-F!"} they can be defined as follows: 
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The same definitions are applied to X"^ and X" , and to all these time series 
for the southern hemisphere. 

3. The periodicities during high- and low-activity periods 

To find the periodicities during high-activity periods I apply the auto-correla- 
tion function (c,-). The original, the positive and the negative fluctuations are 
considered separately. The same is for transformations. I divide each of them 
into 12 parts of length L sa 11 years. For example, for the northern hemisphere 

I consider {i^"}ii, {^nf=i, {^^"+1^1, {^^1^1, {Fr}ti and {Xr}ti 
for each of 12 cycles. Figs. [2^-c show their functions Cr for cycle 18. For the Cr 
of {F"} (Fig. [5^, top) a global maximum is greater than 2(t at t = 11 rotations, 
for {X"} (bottom) the value cn is less than la. Moreover, for {F"^} the value 
Ci2 is between la and 2a and for {X""'"} it is less than la (Fig.[2)D). A decrease 
of Cr at T e [7, 13] rotations, which is at least la, is in 71% of 24 cases (12 cycles 
in each hemisphere). Because the transformation of {F"} into {X^} stabilizes 
the variance of {FT'} such a decrease is probably created by strong fluctuations 
from the high-activity period. Moreover, the width of the Ct maximum of {Fl } 
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(evaluated for the la level) is between two and five rotations. There are also a few 
cases (for example, for {F/"*"} cycles 13, 14, 18, 19, 20) for which each function Cr 
has two or three peaks (their width for the la level is one rotation). This means, 
that the strong fluctuations are quasi-periodical. Fig. [2j: shows that for {i^"~} 
and {Xj"~} in cycle 18 both cn are greater than 2a. This could indicate that 
the fluctuations from the whole cycle create both peaks. In addition, a decrease 
of Cii (it is at least la) could be created by fluctuations from the high-activity 
period. It is also found in 80% of 24 cases. The longer Cr periods are treated as 
subharmonics of about 10-rotation period (IGet ko', '2009p. 

I also apply the Morlet wavelet ( [Torrence and Compo I1998P to six time series 



for each of 24 cases. Fig. [3^ (top) shows the local spectrum of {F^} for cycle 18. 
Black solid contours denote the 95 per cent significance level for detected peaks. 
The global maximum of the integrated spectrum (right) at r = 11 rotations 
is mainly created by fluctuations from the high-activity period. It confirms the 
auto-correlation result (Fig. [1^). Similar wavelet map is for {F[^'^} (bottom), 
but the peak at r = 11 rotations is significant. The global spectrum shows two 
almost the same peaks at t = 11 and 8 rotations. The wavelet map of {X"} (Fig. 
[HJd, top) for the high-activity period is almost the same as for {F[^}. The global 
spectrum for {X""''} is also almost the same as for {F""*"}, but the wavelet 
power at T ~ 11 rotations is significant for {i^"^} while for {Xf^} it is not 
significant. For {i^" } and {X""} in each wavelet map the strongest peak at 
r = 11 rotations is significant (Figs.|3^ and|3jD, middle). It extends in time during 
the high-activity period and dominates the integrated spectrum. Moreover, the 
wavelet map of {X"~} shows this periodicity during the declining portion of 
cycle 18 and smaller power at r sa 8 rotations during the rising portion of cycle 
18. The same was done for all 24 cases. The wavelet and the auto-correlation 
results are similar for r £ [7, 13] rotations (the correlation between them is about 
0.9 for 87%, 72%, 92%, 61% and 75% of the c^ peaks of {F^}, {Ff +}, {Ff"}, 
{X,"} and {Xl'+} respectively). For {Xr} the correlation is about 0.8 for 83% 
of Cr peaks. 

To find the periodicities during the low-activity periods I examine the wavelets 
for 12 cycles together. The map of {X"^} shows the significant peaks at about 
10 rotations from these periods (cycles 13-14, 15-16, 17-18, 20-21). For {Xf"} 
there are also several significant peaks (cycles 13-14, 16-17, 17-18, 18-19 and 20- 
21). Both maps of {X^"^} and {X^^} also show significant peaks during these 
periods (cycles 12-13, 13-14, 15-16, 16-17, 18-19 and cycles 12-13, 15-16, 16-17 
respectively) . 



4. Discussion 

The transformation of fluctuations into the new time series with a constant 
variance enables one to evaluate the time intervals where the functions Cr have 
the most dominant periods. Such analyses indicate that the mean time distance 
between strong fluctuations from solar maxima is about 10 rotations for each 
hemisphere separately. Moreover, the division of {Xi} for each hemisphere into 
two shorter parts: containing data from low-activity periods {{X*}) and from 
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Figure 3. a-c: Wavelet power spectra of a: {F^}, b: {F"~} and c: {F" } mapping a 
time-frequency evolution of about 10-rotation periodicity. Top values of wavelet power are 
denoted by gradual darkening. Black contours denote significance levels of 95 per cent for 
detected peaks. A cone of influence is marked by the dashed region, e-g: Corresponding global 
wavelet power spectra, d: Time series {F"} for cycle 18. h-n: Same as a-g, but {X"}, {X"'~} 
and {Xf+}. 



high-activity periods ({X.f }) enables one to verify possible differences between 
periodicities existing during these periods in 12 solar cycles. 

For each hemisphere the functions c,- of all three these time series have wide 
global maxima at r £ [7, 13] rotations. The functions Ct of {X"} and {X"*} have 
the significant global maxima at r = 11 rotations, but for {X"°} two local max- 
ima at 8 and 10 rotations are between la and 2a. For the southern hemisphere all 
three global maxima are greater than 2a. For {Xf °} this maximum contains two 
almost the same peaks at r = 8, 10 rotations. These results confirm the wavelets 
for all time series. All global spectra have the dominant peaks at r = 10 rotations. 
For {X"^} and {X^^^} these peaks arc much wider and lower than for {X"~} 
and {X^~}. These wide maxima indicate that the time separations between 
fluctuation maxima oscillate around the average length r. This fact could be 
caused by the differential rotation and the different rotation rate of activity 
complexes during their life time (?, ?). In addition, IGetkol (|2004p shows that 
strong positive fluctuations are always created by 2-4 activity complexes which 
have not their largest size at the same time (the maxima of their contribution 
to the monthly Wolf number are shifted from one to three months) . Namely, the 
fluctuations are not strictly periodical, but quasi-periodical. Moreover, because 
the functions Cr (and the wavelets) of positive and of negative fluctuations 
are different, the positive fluctuations, which describe a strong magnetic flux 
emergence, are rather a sequence of pulses following at fixed time intervals than 
harmonic oscillations ( Obridko and Shelting I2007p . In addition, for {F"^} and 
{F^ } the significant wavelet peaks at 10 rotations usually exist during short 
time intervals of high- activity periods and also exist in the maps of {X""*"} and 
{X^~^} during several low-activity periods. 



I also obtained similar results for the monthly Wolf numbers and for the 
group sunspot numbers during high-activity periods (jGetkoi I2006P . For periods 
around 11-12 rotations [Prabhakaran Nayar et aL\ (j2002p showed a significant 
wavelet power of the daily Wolf numbers for short time intervals during several 
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solar maxima. Taking into consideration that the solar activity fluctuations are 
quasi-periodical, the results obtained for daily data and data used in this paper 
are similar. Moreover, an analysis of the auto-correlation of daily sunspot areas 
for solar cycles 10-20 ( |Bogart[ I1982|) leads to the evidence for the existence of 
a few sites of unusually strong solar activity in each solar cycle that persist for 
about 10 rotations and generally rotate with a period of about 27 days. 



5. Conclusions 

The following results have been obtained: 

1. The auto-correlation and the wavelet analyses of the original, the pos- 
itive and the negative fluctuations and their transformations prefer the 
10-rotation quasi-periodicity. 

2. The differences between the auto-corrclation maxima (for r G [7, 13] rota- 
tions) of fluctuations and their transformations suggest that the 10-rotation 
period is created by strong fluctuations from high- activity periods. 

3. The wavelet analysis confirms auto-correlation results for r e [7, 13] ro- 
tations (the correlation is about 0.9 for 87%, 72%, 92%, 61% and 75% 
of the auto-correlation peaks of {Pp}, {Ff+}, {i^""}, {X,"} and {Zr+} 
respectively, and it is 0.8 for 83% of the auto-correlation peaks of {X^~}). 

4. The wavelet maps of {Xf+}, {Xf"}, {X^+} and {X^'} for all 12 solar 
cycles together show several statistically significant peaks at about 10 rota- 
tions for low-activity periods. This means that during these periods small 
quasi-periodical fiuctuations exist. 
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